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Abstract. We report magnetic susceptibility, specific heat, and neutron scattering 
measurements as a function of applied magnetic field and temperature to 
characterize the S = 1/2 quasi-two-dimensional frustrated magnet piperazinium 
hexachlorodicuprate (PHCC). The experiments reveal four distinct phases. At low 
temperatures and fields the material forms a quantum paramagnet with a 1 meV singlet 
triplet gap and a magnon bandwidth of 1.7 meV. The singlet state involves multiple 
spin pairs some of which have negative ground state bond energies. Increasing the field 
at low temperatures induces three dimensional long range antiferromagnetic order at 
7.5 Tesla through a continuous phase transition that can be described as magnon 
Bose-Einstein condensation. The phase transition to a fully polarized ferromagnetic 
state occurs at 37 Tesla. The ordered antiferromagnetic phase is surrounded by 
a renormalized classical regime. The crossover to this phase from the quantum 
paramagnet is marked by a distinct anomaly in the magnetic susceptibility which 
coincides with closure of the finite temperature singlet-triplet pseudo gap. The phase 
boundary between the quantum paramagnet and the Bose-Einstein condensate features 
a finite temperature minimum at T = 0.2 K, which may be associated with coupling 
to nuclear spin or lattice degrees of freedom close to quantum criticality. 
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1. Introduction 

Molecular magnets form intricate quantum many body systems that are amenable to 
experimental inquiry. Their low energy degrees of freedom are quantum spins associated 
with transition metal ions such as copper (spin-1/2) or nickel (spin-1). The interactions 
between these spins can be described by a Heisenberg Hamiltonian of the form 

Ti- = -X X! Jd$ r ■ Sr+dj (1) 
r,d 

Experiments on these materials can test existing theoretical predictions and also reveal 
entirely new cooperative phenomena[U El G2, HI E] . 

If the generic ground state of transition metal oxides is the Neel antiferromagnet 
(AFM), it is singlet ground state magnetism for molecular magnets. When molecular 
units have an even number of spin-1/2 degrees of freedom they typically form a molecular 
singlet ground state. If inter-molecular interactions are sufficiently weak the result can 
be a quantum paramagnet. This is the state that Landau famously envisioned as a likely 
outcome of antiferromagnetic exchange interactions. [6] There are also indications that 
geometrical frustration on lattices containing multiple triangular plaquettes can suppress 
Neel order and/or produce singlet ground state magnetism as is seen for example in the 
3D system Gd 3 Ga 5 0i 2 (GGG) [HE] and possibly in CuHpCl.[9] 

While magnets with an isolated singlet ground state generally have no zero field 
phase transitions, there are at least two field induced phase transitions that terminate in 
zero temperature quantum critical points. The lowest lying excitations in such systems 
typically have spin S — 1, and thus the lower field critical point, H&, occurs as the S z = 1 
triplet branch is field driven to degeneracy with the singlet groundstate at a certain 
critical wave vector. The system becomes magnetized even in the zero temperature 
limit for fields above H&. As the field is increased further, the system is ultimately 
ferromagnetically polarized at a second QCP characterized by an upper critical field, 
H c2 where the field energy overwhelms the spin-spin interactions. Above this field, a 
second spin-gap should open and increase linearly with applied field above H c2 . 

What happens as a function of field and temperature between these critical points 
depends critically on the inter-molecular connectivity. Three-dimensional systems such 
as Cs 3 Cr 2 Br 9 [Tn] and TlCuCl 3 pH CL1 IS E] exhibit a phase transition to long 
range antiferromagnetic order both as a function of applied field and applied pressure. 
For quasi-one- dimensional connectivity the high field phase should be an extended 
quantum critical phase. However, systems that have been examined thus far, such 
as Cu(N0 3 ) 2 • 2.5H 2 [15], [16], and Cu(L — aspartato)(rI 2 0) 2 p2] exhibit long range 
Neel order indicating that they may not be as one-dimensional as originally envisioned. 
Another complication can be multiple molecules per unit cell which lead to staggered 
fields that break translation symmetry elements and preclude a critical phase. Several 
new spin ladder systems have recently been identified that may enable access to a critical 
phase above . [T8l fT9] 
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Two-dimensional (2D) systems are of particular interest since they are at the 
lower critical dimension for Neel as the parent compounds to high temperature 
superconductors. BaCuSi206 for example features spin-1/2 pairs coupled on a two 
dimensional square lattice. [20J An anomaly in the critical exponents and in the phase 
boundary at low temperatures has been interpreted as the result of a dimensional 
crossover from three to two dimensions close to quantum criticality. [21] In quasi- 
two dimensional SrCu2 (603)2 competing interactions play a central role. Its unique 
geometry frustrates inter-dimer interactions, yielding effectively decoupled dimers [22J. 
While the singlet ground state of the material is known exactly, the three finite- 
magnetization plateaus and the critical phases that surround them are more challenging 
to understand. The purely organic quantum AFM F 2 PNNO forms a quasi-2D network 
of spin dimers with a singlet ground state and a 1 meV gap to magnetic excitations. |23j 
The phase diagram includes a magnetization plateau that extends from hqH = 15 
T to fioH = 25 T as well as regions with continuously varying magnetization that 
may indicate unique 2D quantum critical phases. Neutron scattering is required to 
disentangle these complex and unfamiliar high field phases. Unfortunately few of these 
phases can be accessed with current high field neutron scattering facilities, which are 
limited to fields below 17.5 Tesla. 

The material examined here, piperazinium hexachlorodicuprate (PHCC), is a highly 
two dimensional quantum paramagnet with a lower critical field of 7.5 T that is well 
within reach of current high field neutron scattering and an upper critical field of 37 T 
that can be reached in pulsed field magnetization measurements. The large difference 
between the critical fields indicates strong inter-dimer interactions and a potential 
for novel, strongly correlated two-dimensional physics. PHCC, (Cu4Hi 2 N 2 )Cu2Cl6, 
consists of Cu 2+ spin-^ sites interacting predominantly within the crystallographic a- 
c-planes[24, 25J, as illustrated in figure [H Prior zero-field inelastic neutron scattering 
measurements probed the spin dynamics of the low-temperature quantum paramagnetic 
(QP) phase, and established the strongly 2D nature of magnetism in PHCC. The spin 
gap A = 1.0 meV, and the magnetic excitations show dispersion with an 1.8 meV 
bandwidth in the (hOl) plane. There is less than 0.2 meV dispersion perpendicular to 
that plane, indicating a highly two dimensional spin system. The spin-spin interactions 
in the a-c-plane have the connectivity of a frustrated orthogonal bilayer, as shown in 
figure QJc). Bonds 6 and 3 form two extended 2D lattices linked by Bond 1, which thus 
is J_i_ in the bilayer description. Bonds 2 and 8 provide additional frustrated interlayer 
interactions. 

In the current study, we employed inelastic and elastic neutron scattering as well 
as specific heat and magnetic susceptibility measurements to characterize the low- 
temperature phases of PHCC. Unlike 3D coupled-dimer systems, such as TlCuCl 3 , we 
find a significant separation between the phase boundary to long range order (LRO) 
and the crossover line delineating the high-field boundary of the QP regime. This 
allows access to a disordered renormalized classical (RC-2D) regime near the lower 
QCP. Specific heat measurements in the disordered low-field QP phase agree well with 




Figure 1. Crystal structure of piperazinium hexachlorodicuprate (PHCC), 
(C4H 12 N2)(Cu 2 Cl 6 ). (a) View along the a axis, showing well-separated Cu-Cl planes, 
(b) View of a single Cu-Cl plane. Dotted lines between CI sites show potential halide- 
halide contacts, (c) Second view of a single Cu-Cl plane, showing only the interacting 
Cu 2+ spin sites (solid circles). Significant antiferromagnetic spin-spin interactions 
are depicted as light blue (frustrated) and black (satisfied) lines with the same bond 
numbering used in |24j . Line thickness is proportional to the product of exchange 
constant and spin-spin correlation function, | Jd (S • Sd) | . Vectors represent the ordered 
spin structure determined by elastic neutron scattering at T — 1.65 K and /io-ff = 13.7 
T. Dashed lines are the a — c plane unit cell with room temperature lattice constants 
a = 7.984(4) A, b = 7.054(4) A, c = 6.104(3) A, and a = 111.23(8)°, (3 = 99.95(9)°, 
and 7 = 81.26(7)°, space group PI, from 26J. 

series expansions from a quantum bilayer model. Combining elastic neutron scattering 
and field dependent specific heat results allows a determination of the critical exponents 
associated with the transition to LRO. We have previously reported a reentrant phase 
transition at low temperature near H c i wherein the magnetic LRO can melt into the QP 
phase as the temperature is either increased or decreased [27]. We have proposed that 
coupling between other degrees of freedom is important near the QCP, and thus can 
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have a non-negligible effect on the phase diagram in this regime. We also have recently 
examined the LRO transition and found good agreement with a representation of the 
phase boundary away from the reentrant phase transition in terms of magnon Bose- 
Einstein condensation. In this communication, we provide further experimental and 
theoretical information regarding these phenomena, including field dependent neutron 
scattering, magnetization, and specific heat data, and a more detailed theoretical 
discussion of magnon Bose-Einstein condensation in PHCC 

2. Experimental techniques 

Single crystal samples of deuterated and hydrogenous PHCC were prepared via 
previously described methods [21]. Magnetic susceptibility measurements were 
performed at the National High Magnetic Field Laboratory pulsed field facility at the 
Los Alamos National Laboratory using a compensated coil susceptometer. The sample 
environment consisted of a pumped He 3 cryostat located in the bore of a 50 Tesla pulsed 
field resistive magnet. Measurements were performed with the magnetic field along the 
b-, a*- and c- axes on hydrogenous samples with masses between m = 0.7 mg and 
m = 1.4 mg. An empty coil measurement was performed for each individual pulsed field 
measurement of the sample as a direct measure of the sample-independent background. 
Specific heat measurements were performed as a function of magnetic field, fi H < 14 
T, and temperature, T > 1.8 K, on a deuterated single crystal sample, m = 6.76 mg, 
with H || a* using a commercial calorimeter (Quantum Design PPMS). 

Neutron scattering measurements were performed using the FLEX triple axis 
spectrometer at the cold neutron facility of the Hahn-Meitner Institut in Berlin, 
Germany. The sample was mounted on an aluminum holder and attached to the 
cold finger of a dilution refrigerator within a split coil fioH = 14.2 T superconducting 
magnet. A magnetic field in excess of the superconducting critical field for aluminum 
was always present to ensure good thermal contact to the sample. The sample consisted 
of two deuterated single crystals with a total mass of m = 1.75 g coaligned within 
0.5°. The sample was oriented in the (hOl) scattering plane with H || b. Measurements 
were performed using a vertically focusing pyrolytic graphite (PG(002)) monochromator 
followed by a 60' horizontal collimator. Inelastic measurements used a liquid nitrogen 
cooled Be filter before the sample and a fixed final energy of Ef = 3 meV. A horizontally 
focusing PG(002) analyzer was used for inelastic measurements. Elastic measurements 
employed a flat PG(002) analyzer with 60' collimation before and after the analyzer. 
Measurements characterizing the temperature and magnetic field dependence of Bragg 
peaks used a liquid nitrogen cooled Be filter before the sample and neutron energy 
Ei = Ef = 2.5 meV. Measurements of magnetic and nuclear Bragg peaks for determining 
the ordered magnetic structure were performed using the same collimation with a room 
temperature PG filter before the sample and Ei = Ef = 14.7 meV with a vertically 
focusing PG(004) monochromator and a flat PG(002) analyzer. 
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Figure 2. Magnetic susceptibility \ (left axis) and normalized magnetization M/M sa t 
(right axis) (T = 0.46 K only) of PHCC vs. magnetic field for H || b at different 
temperatures. Markers are plotted for every thirtieth data point. The green line in 
the vicinity of H c \ is a fit to M(H), as described in the text. 

3. Experimental results 

3.1. Susceptibility and magnetization 

Pulsed field measurements of the magnetic susceptibility x{H) up to fi H = 50 T with 
H || b are shown in figure El For T < 4 K, the susceptibility has two sharp features. 
The lower field feature marks closure of the spin gap while the upper feature indicates 
the field where M(H) saturates, and the system enters the fully polarized (FP) regime. 
M(H) at T = 0.46 K, obtained by integrating the corresponding x(H) data, is shown as 
a solid line in figure [3 At low temperatures, these features in x{H) indicate the location 
of the critical fields HqH c i ps 7.5 T and HqH c2 ~ 37 T. We note that although there is a 
local minimum in x{H) at fi H = 23.9 T near H m = \{H C \ + H c2 ), and hence a point 
of inflection in the magnetization, there is no indication of fractional magnetization 
plateaus for T > 0.46 K. 

As T increases, the two features in x{H) broaden, and above T = 9 K, they merge 
into a single rounded maximum that moves to lower fields as the temperature is further 
increased. For T > 20 K, x{H) decreases monotonically with H . Despite the thermal 
broadening above T = 5 K, several regimes may be identified from x{H,T). These are 
most clearly seen in figure (3j which shows x{H, T) at fixed field. The rounded peak 
and exponential decrease in x(T) for low magnetic fields is characteristic of a gapped 
antiferromagnet and is similar to prior results obtained at H = for PHCC [24"ll28| 26J. 
As the field increases, the peak in x{T) shifts to lower T, and the low-temperature 
activated behavior vanishes with the spin-gap. As shown for fi H = 10 T in figure 
E(a), for fields fioH c i < [j,qH < 12 T the rounded maximum that signals the onset of 
2D antiferromagnetic spin correlations leading to the gapped state is still visible (near 
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Figure 3. Magnetic susceptibility of PHCC as a function of temperature for H || b. 
(a) liqH < 20 T, (b) 23 T < ^i H < ^i H c2 , (c) H > H c2 . Data were obtained 
by combining all individual pulsed field measurements and averaging over the range 
fio(AH) — 0.5 T for each magnetic field shown. Lines connecting data points are 
guides to the eye. Curves in (c) for T < 8 K are fits described in the text that yield 
the held dependence of the gap above H c2 shown in the inset to (c) . The solid line in 
the inset to (c) has the slope g^sH . 

T = 12 K for fioH = 10 T), but there is then a sharp rise in x{T) as the system enters 
the strongly magnetizable regime at lower T. As shown in figures [3]^a) and [3](b) the 
boundary of this regime continues to move to higher T for fi H < 25 T, before moving 
back to lower T as H approaches H c2 . In the FP region above H c2 activated behavior 
returns, as shown for fi Q H = 40 T and fi H = 47 T in figure M^c). This indicates 
a second spin-gap regime. Note that as the field increases above H C 2, the onset of the 
activated susceptibility moves to higher T, indicating that this second spin gap increases 
with field. An overview of x{H, T) is given in the color contour plot in figure UlTa) and 
will be discussed in more detail below. 

Low-temperature pulsed field susceptibility measurements were also performed at 
T = 0.46 K for H || c and H || a*. These are shown together with the H || b data 
near H c \ and H c2 in figure HI Although the overall shape of x{H) is similar in all 
three orientations, the differences in the onset fields shown in figure H] indicate that 
the critical fields are orientation-dependent. Because the ratio of the critical fields is 
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Figure 4. Magnetic susceptibility of PHCC at T = 0.46 K for H || b, H || c, and 
H || a* in the vicinity of (a) H c i and (b) H C 2- Data obtained for different sample 
masses were normalized to the saturation magnetization for H || b. Solid red lines are 
fits in the vicinity of the upper critical fields as described in the text. Markers are 
plotted every seventh datapoint. 

nearly the same for each of the three single crystal orientations, [H c i(b) / H c2 (h) = 
0.203(3), H cl (c)/H c2 {c) = 0.203(3), and H cl (a*) / H c2 (&*) = 0.205(3)], differences in the 
magnitude of the critical fields are likely due to ^-factor anisotropy associated with the 
coordination of the Cu 2+ ions in PHCC. 

3.2. Specific heat 

Measurements of the specific heat as a function of temperature and magnetic field are 
shown in figure [H Prior measurements were performed on protonated samples in the 
same orientation, H || a*, but were restricted to T < 2 K and ^qH < 9 T. The 
total specific heat C p is shown for representative fields in figure Eta). In zero field, 
C p shows a broad rounded maximum at T pa 7 K, with an exponentially activated 
region at lower T associated with the formation of the spin-gap. Above H c i the form of 
the high-temperature specific heat changes, and as shown for fioH = 14 T, a lambda- 
type anomaly can be seen that is associated with the phase transition to magnetic 
LRO at T = T/v. In figure Mh) we plot the magnetic portion of the specific heat, 
C m = C p — Ci a ttice, where Ci at ticc T 3 is the lattice/phonon contribution to the specific 
heat determined by comparison of the data to the bilayer model discussed later in the 
text. This lattice contribution is shown as a dashed line in figure Eta). Figure [5](b) 
shows the changes in the low-field activated behavior of C m as H increases and the 
spin-gap closes, as well as the increase in magnitude of the lambda peaks as H increases 
above H N . The area under the peak is greater for larger applied fields indicating that 
the T = ordered moment increases with magnetic field. The field dependence of C p 
is plotted in figure [5](c) for selected temperatures. We observe a sharp peak at the 
transition to LRO at H = H^, which has the same field and temperature dependence 
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Figure 5. Specific heat as a function of temperature and magnetic field for deuterated 
PHCC with H || a*, (a) Total specific heat C p /T vs. T 2 at fixed H. Every tenth data 
point is plotted for fioH = 14 T. Calculated lattice contribution is shown as a dashed 
line, (b) Magnetic specific heat C m — C p — Ci at ticc vs. 1/T. Solid lines for fj, H < 10 
T in (a) and (b) are the result of a global fit to the frustrated bilayer model discussed 
in the text, (c) Isothermal specific heat. 



as that seen in C P (T) at fixed H. 
3.3. Elastic neutron scattering 

Neutron diffraction was used to explore the nature of the LRO transition indicated by 
heat capacity data. Field and temperature dependent antiferromagnetic Bragg peaks 
were observed at wave vector transfer Q = ( ( 2n + 1 ) ; 0, ( 2 ™ +1 ) ) f or integer n and m 
throughout the (hOl) plane, corresponding to a magnetic structure with a doubling of 
the unit cell along the a- and c-axes. 

Figure [6] depicts a wave-vector scan through the (|0^) Bragg peak as a function of 
(a) magnetic field and (b) temperature. The scattering intensity at T = 1.7 K increases 
as a function of applied field. At [i H = 14.2 T (figure E](b)) there is only a flat 
background at T = 10.1 K; however, for T = 3.7 K the broad maximum indicates short 
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Figure 6. Wave- vector dependent clastic scattering from PHCC as measured through 
the (|0|) Bragg peak along the (kOl) direction, (a) Magnetic field, H || b, dependence 
at T = 1.7 K. (b) Temperature dependence at fi H = 14.2 T. 




|i H (T) 



Figure 7. Magnetic field, H |j b, dependence of elastic scattering intensity of the 
(ifji) Bragg peak at select temperatures. Temperatures presented increase according 
to the figure legend from left to right. Solid lines are fits described in the text over 
the region /i H < ^H N + 0.75 T. Inset: magnetic order parameter critical exponent, 
(3, determined by fitting a power law 1(H) oc (H — Hn) 2 ^ to magnetic field and 
temperature dependent scattering from the (^0^) antiferromagnetic Bragg peak. The 
Data point at T = 3.69 K corresponds to the temperature scan shown in figure [H 



range AFM correlations. At T = 1.7 K the resolution limited peak indicates long range 
antiferromagnetic order. Scans along the (/i0|) and (C0C) directions were also resolution 
limited. Magnetic Bragg peaks were found at locations that were indistinguishable from 
the nuclear half- wavelength, A/2, Bragg peaks, which were observed upon removal of 
the Be filter to admit A/2 neutrons from the PG(004) monochromator Bragg reflection. 
No indication of incommensurate order was found along the and (|0Z) directions 

through the entire Brillouin zone at T = 1.9 K and [iqH = 11.8 T. 

To characterize the onset of antiferromagnetic order, the (|0|) Bragg peak was 
measured as a function of magnetic field at sixteen temperatures for 0.065 < T < 2.85 
K. Data were obtained by continuously sweeping the magnetic field (podH/dt = 0.01 or 
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Figure 8. Temperature dependence of (|0i) antiferromagnetic Bragg peak elastic 
scattering intensity for PHCC at {IqH = 14.2 T with H || b. Data points were obtained 
from the integrated intensities of rocking curves measured at each temperature. Inset: 
Log-Log plot of integrated intensity as a function of the difference between the 
determined Neel temperature, Tm = 3.705(3) K, and the sample temperature. Solid 
lines are fits as described in text over the range T > T/v — 0.33 K. 



0.03 T/min.) while simultaneously counting diffracted neutrons over a fixed time interval 
of 0.5 minutes per point. At T=0.06 K the measured field trace was indistinguishable 
for the slower and faster rates. The results are shown for select temperatures in figure 
The temperature dependence of the (|0|) Bragg peak at hqH = 14.2 T was measured 
by integrating over a series of rocking curves at 37 temperatures as shown in figure [HJ 
The intensity increases continuously upon cooling below a distinct critical temperature 
as expected for a second order phase transition. 

Figures [9](a) and [9](c) show low-temperature magnetic field sweeps of the (|0|) 
magnetic Bragg intensity plotted as a contour plot and as intensity versus field in the 
vicinity of the critical field for the transition to Neel order, which we denote by H^. The 
data indicate that the ordered phase is thermally reentrant for hqH « 7.5 Tesla. Upon 
heating from 0.065 K, initially decreases with increasing temperature. However, 
there is an apparent minimum in the critical field curve at T = 0.20 K, above which 
temperature Hn increases on heating. The same data are plotted as a function of 
temperature at fixed fields in figure EJ^b). These data show a peak in the scattering 
intensity at T = 0.20 K for fields between fi H ph 7.45 T and HqH ph 7.9 T. The peak 
location remains at T ~ 0.2 K until vanishing for fields above « 8.1 T. 

3.4- Inelastic neutron scattering 

Inelastic neutron scattering was used to probe spin-dynamics as a function of magnetic 
field and temperature. In particular, constant wave vector scans at Q = were 
performed as a function of magnetic field at T = 0.06 K and T = 4.5 K. A portion 
of these data is shown in figure [10J The low-temperature scan at H — T features 
a resolution limited excitation at fno ~ 1 meV, consistent with the previous zero field 
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Figure 9. Reentrant behavior of low-temperature LRO phase as a function of 
magnetic field, H || b, and temperature, (a) Elastic scattering intensity of the (|0|) 
antifcrromagnetic Bragg peak as a function of applied magnetic field and temperature. 
False color image obtained by combining all individual magnetic field order parameter 
scans, (b) Elastic scattering intensity of the (5OI) antiferromagnetic Bragg peak as 
a function of temperature. Data were obtained by combining all of the individual 
magnetic field order parameter scans and averaging over I-IqAH = 0.1 T for each 
magnetic field shown. The symbols on the vertical axis in (a) correspond to constant 
temperature data shown in (b). (c) Magnetic field order parameter scans of the (i0|) 
antifcrromagnetic Bragg peak performed in the region of the reentrant phase transition. 
Solid lines are fits described in the text. 

measurements [24]. In an applied field this excitation splits into three components 
indicating a zero field triplet 1 meV above a singlet ground state. The intensity 
distribution is consistent with the upper and lower peaks resulting from transitions 
to S z = ±1 states while the central highest intensity peak corresponds to a transition 
to an S z = state. The lower energy, S z = +1, peak is observable until « 6 T, 
where it is obscured by elastic nuclear incoherent scattering. Both remaining gapped 
excitations increase in energy continuously as the field is increased above H c ±. The 
higher energy excitation increases in energy and decreases in intensity at a faster rate 
than the excitation which evolves out of the S z — component of the H < H c \ triplet. 
There are also subtle differences in the temperature and field dependence of the width 
of the excitations, which we shall return to later. 
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Figure 10. (a) Constant wave-vector (^Oi) inelastic neutron scattering from PHCC 
at T = 0.06 K (open symbols) and T = 4.5 K (solid symbols) for different values 
of applied magnetic field, H || b. Solid(dashcd) lines are fits to Gaussians for the 
low(high) temperature data as described in the text. Data are offset vertically for 
clarity, (b) Peak locations derived from data such as that shown in (a) as a function 
of applied magnetic field. Solid lines are a fit to the T = 0.06 K data as described in 
text. Dashed lines extend these lines to the origin. 

4. Discussion 

4-1- Susceptibility and magnetization 

A second view of the magnetic susceptibility data shown in figures |2] and [3] is provided 
by the color contour map in figure [TIT a). Since the sharp features in x{H) at the 
lowest-T herald the onset of finite magnetization at H cl and the entry into the fully 
polarized (FP) regime at H C 2, we identify the continuation of these features for T < 5 
K as experimental signatures of the closing of the low-field spin gap and the opening of 
the high-field spin gap with increasing field in the vicinity of the critical fields H c \ and 
H c2 , respectively. These features provide a starting point in constructing an H-T phase 
diagram for PHCC, and are shown as open circles in figure ITIT b). By extending these 
points through inclusion of the sharp features visible in x(T) for T > 5 K, they are 
seen to demarcate the entire large-susceptibility region in figure ITIT a). At non-zero T, 
these points delineate a crossover, rather than a phase transition. Near the lower QCP 
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Figure 11. (a) Color contour map of x(H,T) for PHCC. (b) Magnetic field versus 
temperature phase diagram of PHCC. H c \ and H C 2 are zero-temperature QCPs 
described in the text. Data points were obtained from pulsed field susceptibility 
(x(T,H)), inelastic neutron scattering (INS), elastic neutron scattering (ENS), and 
specific heat measurements (Cp(T, H)). Specific heat measurements were performed 
with H || a*; all other measurements were made with H || b. Values from Cp(T,H) 
have been renormalized as discussed in the text. LRO indicates a long range 
ordered phase; QP is a quasi-2D quantum paramagnetic phase; FP corresponds to 
a ferromagnetic polarized phase; RC-2D indicates a quasi-2D renormalized classical 
phase; and P corresponds to the semi-classical paramagnetic region. Dashed lines 
are guides to the eye. Open diamonds represent the rounded peak locations in x(^) 
delineating a higher temperature classical paramagnetic regime from the correlated 
low-temperature regimes, QP and FP. Open circles represent locations of sharp features 
in x(T) and x(-ff). Closed diamonds correspond to the closing of the spin-gap as 
determined from INS data in figure [T5l Closed squares represent transitions to long 
range order as determined from the field and temperature dependence of the (|0|) 
antiferromagnetic Bragg peak shown in figures [7] and [8] (c) Log-log plot of the low- 
temperature reentrant behavior of the LRO phase. 



H c i, for example, this boundary marks the change from the non- magnetized quantum 
paramagnet (spin-gap) regime to a quasi-2D renormalized classical (RC-2D) region with 
finite magnetization. 
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We also identify two additional crossovers in the H-T plane of PHCC that are 
defined by rounded maxima in x{T) for H < if cl and H > H c2 . These are demarcated 
by open diamonds in figure ITTlfb). At low field, the resulting boundary indicates where 
antiferromagnetic correlations begin to develop as the system is cooled toward the spin 
gap regime. This process has been well-studied in other spin gap systems, such as 
Haldane spin-1 chains [29]. The intersection of this crossover with the zero-field axis in 
figure ITTlfb) marks the onset of a high-temperature semi-classical paramagnetic regime 
(P) where the system has not yet attained significant antiferromagnetic correlations and 
the magnetic susceptibility is that of a disordered ensemble of exchange-coupled spins 
with one energy scale, i.e. Curie- Weiss behavior. 

The high field crossover regime (open diamonds above H = H c2 in figure ITTlfb)) 
is not as well characterized experimentally, but should mark the onset of appreciable 
ferromagnetic correlations as the system is cooled at high field toward the FP state. 
The excitations in the FP state are expected to be ferromagnetic spin waves. As for a 
conventional ferromagnet in a field, these excitations will show a field-dependent gap, 
although here the energy scale for the gap should be set by A FP = g\i B {H — H c2 ). Fits 
of the form x(T) cx e _A/,fesT , such as those shown in figure [3]^c), yield estimates for the 
gap in the FP regime. Values of App as determined from a simultaneous fit of x(T) in 
the temperature range 0.75 < T < 6 K for 40 < /j,qH < 47 T are shown in the inset of 
figure 13(c) and are consistent with a linear increase of App with field above H c2 . 

The evolution of the magnetic susceptibility as a function of magnetic field 
and temperature has been examined numerically for the two-leg spin-ladder ]30j. 
Calculations depict a low-temperature field dependent magnetization, which is 
symmetric about H m = \(H C \ + H c2 ) with square root singularities near the upper 
and lower critical fields, consistent with results for other ID spin-gap systems [31 J . 
Its phase diagram is thus similar to our findings for PHCC. The calculations also 
report similar susceptibility anomalies: a double peaked function at low temperatures, 
a single rounded peak that moves to lower fields as the temperature is increased, and 
a monotonically decreasing field dependence as the temperature is increased beyond a 
certain temperature, T c /. A crossover temperature marking the boundary between the 
classical and quantum regimes, was identified as T c i = 1.109A/&B for the two-leg spin 
ladder. For PHCC, this corresponds to T c i = 1.109A/&B = 13.1 K which agrees well 
with the zero- field lower temperature bound of the crossover regime defined in figure [TTI 
where we find T c i « 15.4 K based upon susceptibility measurements. 

Close to the critical fields the low-temperature field dependent magnetization should 
follow power laws of the form M oc \H — He] 1 / 5 , where H c is a critical field [32]. In terms 
of magnetic susceptibility, this power law becomes 



immediately above (x+) an d below (x-) the respective lower and upper critical fields. 
For PHCC Figures El and |4] indicate that 5 2 > 8\ ~ 1 to account for the observed 



X+ = dM/dH oc (H - H cl ) 
x _ = dM/dH oc (H c2 - H) 



(2a) 
(2b) 
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curvature near the critical fields. Because Si ~ 1, we fit M(H) rather than x{H) to 
determine S. This is done simultaneously for all T = 0.46 K single crystal orientations 
with a common value of Si over a range of one Tesla above H c i. The single crystal 
low-temperature, T = 0.46 K data were also simultaneously fit over a range of one 
Tesla below if c2 using Eqn. (l2~fcj) . The simultaneous fits shared a common S parameter; 
whereas critical fields and overall amplitudes were allowed to vary between the different 
orientations. The resulting values, Si = 0.91(1) and 5 2 = 1-8(2), and their corresponding 
power-laws are shown in figures [2] and HI For comparison, ID singlet ground state 
systems are predicted to have symmetric susceptibility and magnetization curves with 
Si = S2 = 2 [32], and the 2D square lattice Heisenberg model has an analogous upper 
critical field with predictions ranging from <5 2 = 1 to S2 = 1 with logarithmic corrections 
[33~1 [34] . An interpretation of H c i as magnon Bose-Einstein condensation yields the 
prediction of Si = 1 (see section l4~5l) . 

4-2. Specific heat 

We compare our specific heat measurements to the frustrated bilayer model of Gu and 
Shen [35]. Their model contains a strong interlayer bond Jo, a weaker intralayer bond 
Ji, and a frustrating interlayer bond J 2 . This exchange topology is similar to that in 
PHCC. We include a magnetic field dependent term in their thermodynamic expansion 
via a linear Zeeman splitting of the triplet excited state as a function of applied magnetic 
field. We carried out a global fit of all our specific heat data with fioH < 9 T and T < 30 
K using only the free energy terms up to powers of 1/ E~ 2 in Gu and Shen's expansion, 
one independent prefactor for each magnetic field, and global terms consisting of the 
three exchange constants, a ^-factor and a lattice contribution oc T 3 . The results of this 
fit are shown as solid lines in figure [5] and account very well for the temperature and 
field dependence with J = 2.08(1) meV, J x = 0.34(3) meV and J 2 = -0.13(3) meV 

The exchange parameters for the frustrated bilayer model determined from the 
specific heat data for PHCC can be used to evaluate the corresponding dispersion 
relation, equation 4 of Ref . [35] . Figure [12] depicts the calculated dispersion for PHCC 
along several high-symmetry directions in reciprocal space, compared to the zero-field 
triplon dispersion measured by inelastic neutron scattering [2H[25]. The values derived 
from the fits to C m for the spin gap A = 1.11(1) meV and bandwidth of 1.92(2) meV 
agree quite well with the measured spin gap, A = 0.99 meV and bandwidth, 1.73 
meV. Furthermore, the frustrated bilayer model reproduces the shape of the dispersion 
remarkably well given that specific heat data provided the only experimental input. 

The lambda-type peaks in the magnetic heat capacity yield an additional measure 
of the boundary on the magnetic field vs. temperature phase diagram for PHCC, 
figure fTlTb). between the quasi-2D RC-2D disordered phase and the LRO phase. 
However, C p was measured with H || a*, while all the other data included in figure [TlTb) 
were measured with H || b. As noted above, differences in the critical fields measured 
for x{H) for different field orientations are accounted for by g-factor anisotropy. After 
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Figure 12. Zero-field magnon (or triplon) dispersion in PHCC. Data taken from 
Refs [23] and [5S] . Solid lines are the calculated dispersion based upon fitting magnetic 
heat capacity to the bilayer model as described in the text. 



rescaling by the ratio if c i(H || b)/i? c x(H || a*) = 0.917 determined from x{H), the 
phase boundary derived from C p (H, T) is seen to agree very well with that determined 
from elastic neutron scattering as shown in figure fTTIfb) and (c) and discussed below. 

The magnetic field dependent specific heat shown in figure EJc) also provides a 
sensitive measure of critical scaling in the vicinity of the LRO transition. C m (H) 
is shown in figure [13] as a function of reduced magnetic field, h = | —jf^ | near the 
transition field H^. Power law behavior is clearly seen for H > H^, while below H^, 
the magnetic specific heat asymptotically approaches the same power law behavior as 
the H > H N portion of the data. In addition, there is a region of reduced field for 
H < Hn, 0.02 < h < 0.2, where we observe a different power-law behavior, i.e. a 
different slope on the log-log scale plot as indicated by the dashed line in figure [13] The 
two different power-laws observed for H < H^, as shown by the fits for 0.02 < h < 0.2 
and the asymptotic change in power-law as one approaches Hn is a potential indication 
of an additional disordered portion of the phase diagram for H < H^. This is in fact 
directly measured as a gapless disordered region through comparisons of the magnetic 
susceptibility and neutron scattering measurements as we discuss later in the text and 
as shown in figures [TlTb) and (c). 

We fit the H > data in figure [13] to a power law of the form 

C m = A(h)~ a , (3) 

where a is the specific heat exponent. A global fit for T = 2 K and 3 K results in 
a = 0.097(2) while individual fits yield a{T = 2K) = 0.104(3) and a(T = 3K) = 
0.072(2). Individual fits to Eq. [3] for H < H N over the range 0.02 < h < 0.2 yield 
a(T = 2K) = 0.175(2) and a(T = 3K) = 0.2626(7). These values are larger than those 
obtained for H > H^, which indicates that the fits include data beyond the critical 
regime. Table [1] summarizes theoretical and experimental results for the 3D magnets. 
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Figure 13. Magnetic heat capacity as a function of reduced magnetic field, h, in the 
vicinity of the phase transition to LRO. H < data (solid symbols) are plotted 
against h = H x~ H ■ H > Hn data (open symbols) are plotted against h — H h^ n ■ 
Black solid lines correspond to a global fit of H > Hn data to a power law as discussed 
in text. Coloured solid lines correspond to individual power law fits for H > Hjy. 
Individual fit results are very similar to the global fit result. Dashed lines correspond 
to individual power law fits for H < Hn for 0.025 < h < 0.2 in the RC-2D regime. 

Table 1. Measured exponents for PHCC compared to calculated exponents listed 
in [Ml 133 HO Hi] for both 3D Heisenberg (HB), XY (3D), and Ising (2D and 3D) 
transitions. Errors in referenced values are listed when available. Error bars for a 
in PHCC correspond to one standard deviation greater and less than the respective 
maximum and minimum fitted values described in the text. 

Calc. HB Calc. XY 2D Ising Calc. 3D Ising 



PHCC Ref |3Z] Ref [35] Ref [36] Ref [33] 



a 


0.09718:8? 


-0.1336(15) 


-0.0146(8) 





0.1099(7) 


P 


0.34(2) 


0.3689(3)* 


0.3485(2)* 


0.125 


0.32648(18)* 


7 


1.22(4)* 


1.3960(9) 


1.3177(5) 


1.75 


1.2372(4) 


<5 


4.6(2)* 


4.783(3)* 


4.780(2)* 


15 


4.7893(22)* 


V 


0.634(1)* 


0.7112(5) 


0.67155(27) 


1 


0.63002(23) 


V 


0.07(6)* 


0.0375(5) 


0.0380(4) 


0.25 


0.0364(4) 



* Calculated from scaling relations and other listed critical exponents. Error bars 
propagated from the measurements or calculations via the scaling laws. 

4-3. Elastic neutron scattering 

4-3.1. Ordered magnetic structure To determine the field induced magnetic structure 
we measured fourteen magnetic Bragg peaks in the (hOl) plane. The measurements 
were carried out at T = 1.65 K and fioH = 13.7 T which is within the ordered phase 
as established by heat capacity measurements. Table [2] shows raw experimental results 
in the form of wave vector integrated intensities, J m (r m ), for transverse (rocking) scans. 



Field-driven phase transitions in a quasi 2d quantum AFM 19 
I m (r m ) is related to the magnetic structure of the material as follows: 

I m {r m ) = C / ft T (q ± )dq ± iV m ^^J m (r m ). (4) 

Here C is a normalization factor that depends on aspects of the instrument such as the 
incident beam flux that are not considered in the elastic resolution function 7£ T (q)[40j 
which is sharply peaked around q = 0. N m and v m are the number of magnetic unit 
cells and the magnetic unit cell volume respectively. The magnetic Bragg scattering 
intensity is given by 

J m (Q)=r 2 (|^(Q)| 2 -|Q-^(Q)| 2 ), (5) 

where ro = 5.38 fm.j^T] 

^(Q) = ET f d(Q))e- H ' d(Q) S d (6) 
d z 

is the magnetic structure factor, g& is the Lande g-factor for the magnetic ion on site 
d within the magnetic unit cell, and i*a(Q) is the magnetic form factor of the Cu 2+ 
ion (42J The magnetic structure enters through the expectation values for the spin 
vector, Sd in units of h. 

Table [2] compares the measured scattering intensity I m {j m ) to the calculated 
intensity for the uniaxial magnetic structure that best fits the data and is shown in 
figure [U The ordered moment lies within ~ 7(4)° of the a-c-plane and the in-plane 
component is within ~ 6(4)° of the a-axis. For T = 1.65 K and fi^H = 13.7 T 
the sub-lattice magnetization is gfisiS) = 0.33(3)g/ie, determined by normalization 
to incoherent scattering of the sample. Assuming antiferromagnetic interactions, the 
ordered structure is consistent with the sign of the bond energies determined on the 
basis of zero-field inelastic neutron scattering. [21] 



4-3.2. Phase diagram Figures [3, [8] and [9](c) compare the temperature and field 
dependence of the magnetic Bragg intensity at Q — (|0|) to power laws 1(H) oc 
(H—H^) 213 and I(T) oc (T^—T) 213 in order to extract the critical exponent, (3. Rounding 
in the vicinity of the phase transition was accounted for as a distribution in the critical 
field. Possible origins of the 3.6% distribution width include inhomogeneities in the 
applied field and chemical impurities. Least squares fitting included only data within 
0.75 T and 0.33 K of criticality. The corresponding fits appear as solid lines in the 
figures. 

Values of (T/v, Hn) extracted from the analysis appear as filled squares on Figure fTTl 
There is a clear distinction between the onset of a state of high differential susceptibility 
as determined from x{H, T) and the boundaries of the long range ordered phase (LRO) 
extracted from neutron diffraction data. We denote the state of the system between 
these boundaries as 2D renormalized classical. Spin correlations in this phase were 
examined by quasi-elastic neutron scattering measurements along the (|0Z) and (/i0|) 
directions through the (|0|) wave-vector at T = 3 K and fi H = 9 T. No magnetic 
Bragg peaks were encountered. Instead the broad peak centered at Q = (^0|) indicates 
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Table 2. Measured and calculated relative magnetic Bragg peak scattering intensities 
at the indicated magnetic reciprocal lattice points r m = ha* + Ic* . The parentheses 
indicate statistical errors only. The minimum in the crystallographic reliability index 
R = \Im S ~ ^m lc \/J2^m S = 0-3 is found for the magnetic structure depicted in 





[T] and described in 


UllVj V V.' W L ■ 


h 


I 


I m (obs) 


I m (calc) 


1 


-1/2 


0.00(1) 





1 


1/2 


0.00(1) 





-1 


-1/2 


0.00(1) 





1/2 


-1 


0.00(1) 





-1/2 


-1 


0.00(1) 





3/2 


1/2 


0.01(1) 


0.002 


1/2 


-1/2 


1.00(2) 


1 


-1/2 


-1/2 


0.21(1) 


0.234 


3/2 


-1/2 


0.16(1) 


0.333 


1/2 


-3/2 


0.11(1) 


0.250 


5/2 


1/2 


0.03(1) 


0.027 


7/2 


5/2 


0.04(2) 


0.154 



short range correlations with the same local structure as the long range order that sets 
in at slightly higher fields. 

4-3.3. Critical exponents The order parameter critical exponent, /3, determined for 
each of the fits to the magnetic field and temperature dependent scans is shown versus 
temperature in the inset to figure [7J For T > 0.5 K, the average value is /3 — 0.34(2) 
which is indistinguishable from the critical exponent for the 3D Heisenberg [37], Ising 
(P = 0.326), and XY (/3 = 0.345) mo dels [36]. However, the positive specific heat critical 
exponent, a, uniquely identifies the phase transition as being in the 3D Ising universality 
class. This is consistent with expectations for a phase transition to long range magnetic 
order in space group PI. We also calculate the remaining critical exponents based upon 
the Rushbrooke , a + 2/3 + 7 = 2, Fisher, 7 = v{2 — 77), Josephson, 2 — a = vd, and 
Widom, S = 1 + Z, scaling laws, assuming d = 3 dimensions [361 S3]- The values are 
again consistent only with those of 3D Ising transitions as summarized in table [U 
Recalling the exponents 81 = 1.022(2) and 62 = 1-8(1) in (!2al) and (j2~5jl determined 
from pulsed field magnetic susceptibility data, note that these values are different from 
S derived from the scaling relations. This is because 5± and 62 are associated with the 
uniform susceptibility while 5 in table table Q] characterizes the response to a staggered 
field conjugate to the Neel order parameter. 

The apparent order parameter critical exponent increases as temperature is reduced, 
such that (3 = 0.40(1) for T < 0.43 K (see inset to figure[7J). A possible interpretation of 
this behavior is that the system enters a crossover regime associated with the approach 
to quantum criticality. Indeed, in the T = limit the system should be above its upper 
critical dimension such that mean field behavior with a critical exponent (3 m f = 0.5 
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should be expected. 

The field sweeps of the Q = (|0|) magnetic Bragg peak intensity shown in 
figure E(a)-(c), indicate a reentrant phase transition. The H — T phase boundary 
obtained by fitting these data (figure fTlT b)) show a local minimum in the phase boundary 
between the LRO and QP phases. A possible origin of this anomaly is coupling to 
non-magnetic degrees of freedom such as nuclear spin or lattice distortions that are 
unimportant at higher temperatures but become important sufficiently close to the 
quantum critical point [27J. 

4-4- Inelastic neutron scattering 

4-4-1 ■ Gapped Excitations Figure [10] shows the magnetic field and temperature 
dependence of the excitation spectrum at the critical wave vector Q = (|0|) for < 
H < H C 2. We fit each spectrum to a sum of Gaussian peaks to account for the magnetic 
excitations with a constant background and a field independent Gaussian+Lorentzian 
lineshape centered at zero energy to account for incoherent elastic nuclear scattering. 
The field dependence of the peak locations are plotted for the T = 0.06 K and T = 4.5 
K data in figure fTDTb). The data demonstrate the triplet nature of the excited state 
through its Zeeman splitting. Above H c i, the rate of change of the peak position 
as a function of field for the higher energy peak is twice that of the lower energy 
excitation, indicating a magnetized ground state. We fit the low T field dependent 
peak locations to the spectrum of an antiferromagnetically coupled spin-1/2 pair with 
singlet-triplet splitting A subject to a magnetic field described by Hz = —g^B^H- 
Results of this two-parameter fit to the low-temperature peak locations are plotted as 
solid lines in figure fTDT b). and yield values of A = 1.021(2) meV, g = 2.33(1), and 
HoH c i = A/(<7/iB) = 7.58(3) T, which are in agreement with zero-field inelastic neutron 
scattering and bulk magnetization data. 

The average resolution-corrected half width at half maximum, T, of the three 
Gaussians that were fitted to each spectrum for T = 0.065 K and T = 4.5 K are shown 
in figure fTH While T, considering systematic uncertainties in the energy resolution, 
is indistinguishable from zero at low T, it increases with applied field at elevated 
temperature. Also increasing with field and temperature is the magnon density. Hence 
the result is consistent with an increased rate of magnon collisions as their density 
increases. Self consistent mean field theories have been developed to account for these 
effects. |B] 

The inset to figure fTH shows the field dependent ratio of the integrated intensities 
of the central Gaussian peak to that of the satellite peaks. Based upon the polarization 
factor and the Wigner-Eckart theorem as applied to a Zeeman split triplet, this ratio 
should be 0.5 for H < H c \ jl5]. Both the low- and high-temperature constant wave- 
vector scans are consistent with this result. Above H c \, Bose condensation modifies the 
ground state and the experiment shows a decrease of intensity for the highest energy 
mode. 
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Figure 14. Characteristics of constant energy scans as a function of applied magnetic 
field at T = 0.06 K (open symbols) and T = 4.5 K extracted through Gaussian fits. 
Representative scans are shown in figure [TDJ (a) Average resolution corrected half 
width at half maximum for finite energy modes. Inset shows the ratio of integrated 
intensities for the field dependent modes to the central field independent mode. The 
solid line indicates the theoretical ratio of 0.5 expected for H < H c \ at T = 0. 



We use inelastic neutron scattering to examine the crossover between the low-field 
gapped (QP) and gapless (RC-2D) portions of the phase diagram. The T = 0.06 K lower 
critical field determined by comparing the peak locations in low-temperature constant 
energy scans to Eq. ?? is plotted in figure fTlTa) and (b) (solid diamond at T = 0.06 
K). With wave vector transfer at Q = (|0|) and energy transfer fixed at Tiu = 0.3 meV 
the data in figure [15] show the field dependent intensity for various fixed temperatures. 
At the lowest temperature a sharp peak indicates the field at which the S z = 1 mode 
passes through the 0.3 meV energy window of the detection system. With increasing 
temperature and hence increasing magnon density not only does the peak broaden as 
expected from an increased relaxation rate, it also moves to higher fields indicating a 
finite slope for the crossover line in the H — T plane. Just as for the finite magnon 
lifetime, the origin of this effect is magnon repulsion. The increasing background is 
likely due to population of low-energy gapless excitations which develop for magnetic 
fields above the lower critical field. By extrapolation to the Tiu = line in Fig. [ToT b). 
these data were used to determine the temperature dependence of the field required to 
close the spin gap and the corresponding data points are included in figure fllTa) and 
(b) (solid diamonds). The coincidence between the crossover field determined in this 
manner and the fields associated with differential susceptibility anomalies establish a 
not unexpected association between effectively closing the gap at the critical wave vector 
and maximizing the differential uniform susceptibility. 

4-4- 2. Magnetic Excitations in the high field paramagnetic state The previous results 
indicate a well defined crossover in PHCC from a low field and low temperature state 
dominated by a gap in the excitation spectrum to a high field state with an excitation 
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Figure 15. (a) Inelastic neutron scattering intensity at hu> = 0.3 meV and (|0|) 
as a function of applied magnetic field, H || b. Solid lines are the result of a global 
fit to a Gaussian peak with a sloping background above the magnetic field Hq and a 
constant background below the magnetic field Hq. Values of the constant background, 
Hq, and slope were global parameters in the fit. (b) Energy transfer versus magnetic 
field for peak locations as determined from (a) and the T = 4.5 K and T = 0.06 K 
data in figure HOf a) . Symbols used in (b) correspond to those used in (a) . Solid lines 
are fits to lines as described in the text. The extrapolated intercepts define an effective 
temperature dependent crossover field. 



spectrum that has finite spectral weight at zero energy. The crossover clearly occurs 
at lower fields than required to induce long range magnetic order. There are possible 
analogies to certain 5 = 1/2 ID Heisenberg antiferromagnetic systems including the 
linear chain, the next nearest neighbor frustrated chain, and the alternating chain 
where a low-energy continuum of excitations is predicted to exist upon closure of 
the spin-gap 061 07]. Calculations examining both gapped and gapless 5=1/2 
spin-orbital chains describe an incommensurate gapless continuum for large portions 
of exchange and orbital interaction parameter space [H]. A characteristic of these 
continua are soft incommensurate modes moving through the Brillouin zone from 
the antiferromagnetic zone center to the ferromagnetic point as a function of applied 
magnetic field [2H 021 EH EI]- Incommensurate magnetic field dependent soft modes 
are also found in numerical calculations for the 2D antiferromagnetic square lattice with 
next nearest neighbor interactions, a system with similar antiferromagnetic interactions 
as PHCC 

An important distinction between PHCC and the systems mentioned above lies 
in the fact that for PHCC the putative gapless state so far was only observed at 
finite temperatures. Ascribing this to the inevitable effects of inter-plane coupling, 
the analogue of which occurs via interchain coupling in quasi-one-dimensional systems, 
there exists a potential for field induced incommensurate correlations in PHCC. Here 
we present preliminary experiments at the highest fields presently accessible to neutron 
scattering. The field was fioH = 14.2 Tesla and the temperature chosen was T = 
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Figure 16. (a) Magnetic inelastic scattering from PHCC at the critical wave vector 
(|0|) for fioH = 14.2 Tesla and T = 4.5 K which is immediately outside the long 
range ordered phase. A presumed non-magnetic background measured at (^01) has 
been subtracted. The solid line is a fit to a Lorentzian relaxation function as described 
in the text, (b) Constant energy scans, fito = 0.36 meV, along the (5OZ) direction at 
T = 4.5 K and fi H = 14.2 T for Tiu> = 0.36 and 0.55 meV. Solid lines through data 
are fits to a single peaked Lorentzian functions with sloping backgrounds 



4.5 K, which is 20% above the critical temperature at that field (Tc{hqH = 14.2 
T)=3. 705(3) K). Shown in Fig. [16] (a) are difference data indicating excess magnetic 
spectral weight at the antiferromagnetic zone center (|0|) as compared to (|0l). 
The solid line shows a fit to the following resolution convolved Lorentzian relaxation 
spectrum, 

<->(^J OC ; — — r— -, (7) 

which obeys detailed balance and features a relaxation rate, HT = 0.025 ± 0.2 meV. 
The spectrum extends to and is approximately symmetric about zero energy, which 
is consistent with a 2D renormalized classical phase close to the phase transition to 
commensurate magnetic order. 

Figure [16] (b) shows two constant energy scans at hu = 0.36 meV and hu = 
0.55 meV that indicate short range commensurate magnetic correlations. The fitted 
Lorentzian peaks, which are much broader than the experimental resolution (solid bar) 
correspond to a dynamic correlation length scale of 3.2 ± 2.5 A and 4.2 ± 2.1 A for 
the hu = 0.55 meV and hu = 0.36 meV measurement respectively. The magnetization 
of PHCC is only 20% of the saturation magnetization in these measurements so any 
incommensurate shift, which should be a similar fraction of the distance to the Brillouin 
zone center is not clearly visible in the present data. 

4-5. Bose condensation of magnons 

In this section the phase transition to long range magnetic order in PHCC is analyzed 
from the perspective of magnon condensation [52], [53], IS] [55]. Following Nikuni et al. 
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and Misguich and Oshikawa [22] we treat magnons as bosons whose density n is 
proportional to the deviation of the longitudinal magnetization from its plateau value 
M z q: gfiBn = \M Z — M z q\. The low-density magnon gas is described by the Hamiltonian 

H = 53(e k - /i)4a k + a i+k' a q-k' a q-ka q + k , (8) 

k V q,k,k' 

where \x = gfi B {H — H c i) is the boson chemical potential and Vq is a short-range 
repulsion. The energy dispersion ek is assumed to have the global minimum at 
k = [|0|]. For convenience, excitation energies are measured from the bottom of 
the band, so that ek = 0. 

At zero temperature the system acquires a condensate (a ) = ay/V when /i > 0. 
The energy density of the ground state in the presence of a condensate (neglecting 
quantum corrections) is fo|«| 4 — /^M 2 , which is minimized when \a\ 2 = h/2vq- In this 
approximation the boson density is 

n=\a\ 2 = g/i B (H - H c0 )/2v . (9) 

n is related to the longitudinal magnetization 

M z = g^ B n (10) 

At a finite temperature we obtain an effective Hamiltonian 

H = ^2(e k - fi + v n)ala k , (11) 

k 

where n is the density of thermally excited bosons. The bosons condense when the 
renormalized chemical potential fi = fi — 2vqu vanishes. This yields the critical density 

'^>4? expMT)-r < 12 > 

The critical field is then 

H C (T) = H e (0) + 2v n c (T)/ gf i B . (13) 

The boson scattering vq can be extracted from low-temperature magnetization curves. 
The magnon density per unit cell n is equal to the ratio M/M sat (there is exactly one 
magnon per unit cell when M = M sat ). Based on the magnetization curve at the lowest 
available temperature T = 0.46 K we obtain vq = 1.91 meV x the unit-cell volume 
(u.c.v.). In other words, if we model the magnons as bosons on a cubic lattice with 
anisotropic hopping and on-site interactions U, the interaction strength would be U = 2 
meV. The linear behavior extends over a range of some 10 T, or 13 K in temperature 
units. 

Using this value of the boson coupling we can test the predicted relation between 
the critical field and critical magnon density (fTBl at finite temperatures from the data 
for n = M/M sat in the temperature range 0.46 K < T < 6 K. From these we extracted 
pairs (M, H) whose H(T) values were sufficiently close to the line of phase transitions 
H C (T). The resulting graph is shown in Fig. [T71 The zero-temperature critical field 
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H c i = 7.41 T was determined from the onset of Neel order at the lowest temperatures 
(T < 0.3 K). The agreement with Eq. (|T3|) is fair. 

The line H C (T) is obtained from Eqs. (1121) and (1T3"|) . The magnon dispersion is taken 
from Stone et al. [23] (the 6-term version). That dispersion relation e 2 , is strictly two- 
dimensional, which precludes Bose condensation at a finite temperature. We therefore 
added weak tunnelling between the planes: 

e k = e2?+2<r[l-cos(V»)]. (14) 

The experimental data [57] give an upper bound for the total bandwidth in the [010] 
direction: 4r < 0.5 meV. 

When T < r, only the bottom of the magnon band has thermal excitations and 
one may replace the exact band structure with a parabolic dispersion to obtain 

3DtA 3/2 

v 2nh 2 / 

where m 3D = (m a mftm c ) 1//3 is the 3D effective mass. In the limit of a very weak 
interplane tunnelling we can also access a regime where t <^T ^W, where W ~ 2 meV 
is the magnon bandwidth. In this regime the in-plane dispersion can still be treated as 
parabolic. This yields the critical density for a quasi-2D Bose gas [58] 

m 2D T 2T 
2nTi b t 

where m 2D = (m a m c ) 1//2 and b is the interplane distance. In the intermediate regime 
T = 0(t) the critical density first rises faster than T 3 / 2 [56J before crossing over to the 
TlogT behavior. 



n c (T) = C(3/2) ( y , (15) 
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Figure 18. Squares are data for H C (T) from elastic neutron scattering. H c i{T = 0) 
7.41 T, vq = 1.91 meV u.c.v. Solid lines are discussed in the text. 



Curves for n c (T) were obtained by integrating Eq. ffl2|) numerically for several 
values of the interplane tunnelling r. Since Vq is already fixed, there are no further 
adjustable parameters. The data below 2.5 K (H c < 9 T, n c < 0.05 per unit cell) are 
best described by the curves for r = 0.03(1) meV. The results are shown in Fig. [TBI 

The theory based on the Hamiltonian (fTTl) assumes that the boson interaction v o 
is weak. However, it is also applicable to strongly interacting bosons, as long as the 
interaction is short-range and the gas is sufficiently dilute. In that case vq is not the 
bare interaction strength (which can be large) but rather the scattering amplitude (the 
t-matrix) for a pair of low-energy bosons. For hard-core bosons on a cubic lattice 
with a nearest-neighbor hopping amplitude t, one obtains t>o ~ St [59]. In our case 
t « W/8 = 0.25 meV 2 meV. Another estimate along these lines, using the 

"exact" magnon dispersion yields, via 



yields a similar value of vq = 2.5 meV. These numbers demonstrate that the inferred 
interaction strength of 1.91 meV u.c.v. is quite reasonable. 

At higher temperatures the agreement worsens. At T = 3.7 K — the highest 
temperature where the condensation is observed — the mean-field theory predicts H c = 
11 T, well below the observed value of 14.2 T. The observed magnon density at the 
critical field is 0.19 per unit cell (the theoretical estimate is 0.10). It is possible that 
the mean-field theory of magnon condensation, appropriate for a dilute Bose gas of 
magnons, becomes unreliable due to overcrowding. (At a technical level, the assumption 
of pairwise interactions may break down: two particles no longer collide in a vacuum.) 
Another possible reason for this failure is an increased strength of (thermal) critical 




(17) 
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fluctuations. 



4-6. Kosterlitz-Thouless scenario 

If the planes in PHCC were fully decoupled from each other the Bose condensation would 
change into the Kosterlitz-Thouless (KT) vortex-unbinding transition. It does not seem 
likely that the data can be reconciled with this scenario: in d = 2 spatial dimensions 
the crossover exponent <fi = 1, so that the line of KT transitions H C (T) = H c (0) + CT^ 
is expected to be linear. 

The quantum criticality of the two-dimensional Bose gas was recently examined by 
Sachdev and Dunkel [60] who used a combination of scaling arguments and numerics to 
determine H C (T) for a weakly interacting Bose gas pip . (Analyzing the RG flow alone 
does not solve the problem: the KT fixed point lies in a strongly coupled limit.) 

The RG analysis focuses on two variables: the (renormalized) chemical potential ft 
and the two-particle scattering amplitude v. Unlike in d = 3 dimensions, the scattering 
amplitude vanishes in the limit of zero energy. Therefore the role of the coupling constant 
is played by the t-matrix for particles with energies 0(T): 

4nh 2 

V ~ mlog (A/T) ' ( } 

where A is a high-energy cutoff on the order of the magnon bandwidth. 

While both ft and v flow under RG transformations, a certain combination of these 
variables, namely Tv/ft, remains unchanged. The value of this constant determines 
whether the system flows to the disordered phase or to the phase with long-range 
correlations. The KT transition corresponds to 

^ = -K » -34. (19) 

At the KT transition, the renormalized chemical potential is 

ft = - (20) 
K log (A/T) 

We see that \ft\ is a bit below T (particularly if T A), which means that there are 
many thermally excited particles. 

The magnetic field is connected to the bare chemical potential: 

gLL B {H - H c0 ) = fi= ft + 2vn. (21) 

The density is computed for an ideal gas: 

r d 2 k 1 mT 1 mT T 

U ~J (2tt) 2 e ( ft2fc2 /2 m -A)/T _ 1 " g 1 - e»/ T ~ 2tf? g ^ ( } 

This yields the final result for H C (T): 

AT 



gLi B (H - H c0 ) = ji 



K A\ 2tt 

log( ^ log rJ-¥ 



(23) 



log (A/T) 

It depends weakly on the ultraviolet cutoff A. For A = 2.5 meV, we have li ~ 1.5T in 
the temperature range of interest. The curve is labelled 'KT1' in Fig. [TSJ 
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The result (T2"3"j) can also be obtained directly from Eqs. (3) and (4) of Sachdev and 
Dunkel [60J if one takes the classical limit /i C T. The value \x ~ 1.5T suggests that 
we are (just) outside the classical regime and one needs to solve the Sachdev-Dunkel 
equations numerically. The resulting curve (for A = 2.5 meV) is shown in Fig. [18] as 
'KT2'. It is also approximately linear, with a slope dfi/dT w 1. 

Even with the refined approach, the KT theory does not seem to agree with 
the data. As anticipated, it produces a linear dependence H C (T) (up to logarithmic 
corrections) instead of a higher power law. The rise of the data above the KT line at 
high temperatures may signal the inadequacy of the continuum approximation at high 
magnon concentrations. 

5. Conclusions 

Originally thought to be a weakly alternating spin chain, our measurements establish 
PHCC as a quasi-two-dimensional spin-dimer system. The gap is sufficiently small to 
allow access to the lower critical field and the bandwidth sufficiently large to produce 
several distinct phases. 

The long range Neel ordered phase is a commensurate antiferromagnet with 
moments perpendicular to the applied field direction. The critical exponents for the 
phase transition are similar to those of a 3D Ising model. The transition into this 
phase can be described as magnon Bose-Einstein condensation. Mean-field (Hartree- 
Fock) theory does an adequate job describing the magnetization M(H) at the lowest 
temperatures (T = 0.46 K) and the shape of the critical line H C (T) for T < 3 K. 
The theory predicts a small but nonzero magnon dispersion of 0.08 > 4r < 0.16 meV 
perpendicular to the 2D planes. The higher power of the H C (T) dependence is caused 
by deviations of the magnon density of states from the E 1 ^ 2 behavior of the continuum 
theory. The extra DOS (due to a van Hove singularity) is positive, hence an increase 
over the T 3 / 2 law. The continuum KT theory gives an approximately linear dependence 
H C (T) and thus does not seem to do a good job. 

Close to the lower quantum critical point there is an anomaly in the phase boundary 
that is not accounted for by the BEC theory (see figure [TTTb) and figure E^a)-(c).) 
Low T anomalies near the zero temperature termination of phase boundaries are also 
seen experimentally in the quantum liquids 3 He[6T] and 4 He[62], [63], the frustrated 
spin system GGG[7J [8], and in numerical and theoretical results for the BCC Ising 
antiferromagnet [64], the Ising antiferromagnet on a cubic lattice [65J, and 5 = 1/2 
XXZ ID- (66] and 2D- [67] antiferromagnets. In 4 He the anomaly is associated with 
there being more low energy degrees of freedom in the solid than the liquid and a 
low— T anomaly in the phase boundary of LiHoF4 has been associated with nuclear spin 
degrees of freedom. [68j For PHCC we speculate that features beyond a spin model such 
as the relatively soft lattice or nuclear spins become relevant sufficiently close to the 
quantum critical point. Such a mechanism which appears to be rather general merits 
further exploration through experiments that probe such auxiliary degrees of freedom 
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directly. 

The high field paramagnetic phase is perhaps the most intriguing aspect of PHCC. 
The crossover from the quantum paramagnet to this phase remains remarkably sharp 
upon heating and the experiments show that a finite temperature pseudo spin gap closes 
there. Further analysis is required to understand why this crossover remains so distinct. 
A proper examination of the gapless high field phase requires inelastic neutron scattering 
magnetic fields of order 30 Tesla where the system is more heavily magnetized. This 
work must await development of a higher field facility for inelastic neutron scattering. 
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